هی ی هزیر 9 :10۲ سال سی و یکم» شمارة چهار, ۱۳۹۷ 
استنتاج معادلات جدید برای تخمین بیشینه شتاب و سرعت زمین ناشی از زلز له * 
در حشان/ ٩(‏ صایر ی () 
علی درخشانی علی صابری 


چکیده در این پژوهش یک روش نوین موسوم به 15 برای به‌دس تآوردن معادلات پیش‌بینی حرکات نیرومند زمین مورداستفاده قرار 
گرفت. بيشینه شتاب زمین و بیشینه سرعت زمین بااستفاده از بزرگی زمین‌لرزه فاصلةٌ منبع تا سایت. متوسط سرعت موج برشی و مکانیسم 
کسلش و نیز بهره‌گیری از پایگاه داده‌های کسترده پژوهشگاه مهندسی زمین‌لرزه (ل00) فرمول‌بندی شدند. برای اعتبارسنجی, مدل حاصل با 
مدل‌های شناعته شده مقایٍ سه گردید. تحلیل <سا سیت وآنالیز پارامتریک برای تعیین اهمیت پارامترهای موثر بر مدل و <سا سیت آن به 


تغییرات پارامترها انجام شد. معادلات پیشنهادی بسیا رآسان است و می‌توان آنها را با اطمینان برای اهداف پیش‌طراحی استفاده نمود. 


واژه‌های کلیدی حرکات نیرومند زمین. معادلات پیش‌بینی. مدل درختی ۷15 ارزیابی خطر لرزه‌ای. 


0 ۳۵۵۱۰ ۱۸0۵۵60 ماه مه ۵۲ ممت)متاع! ما مممتاهنو۲ ع۱ ۵ ممتاه۱(۵۲۱۲۵ 
۲ 200 «0۵00عع جر 


ت06 .خر تصقطوحلهرع(1 ,خر 


0 ۳0۱۵۱۵-۱۱0107۱ 6۳۳۷۵ 10 ۵۱۴۵0۵0 ک۱۵۵ رد۷ 6۵0 ۱۱۵۵4 ۵۷ ۵ ای کنا ۲۰ ]026وهاظ۸ 
۷۸۵۵۳۰۱۱۵۵۵۵۵ (2۲۷)ظ) ۷۵۱06 2۲۵09۵۱ 2601 0۵0 (4ی)(ظ) 066۵2۵۵ 2۲۵9۵۱۵ 01( (ع۷()) دوه 
-61607 ۷6۵86 ,015101166 516 10 50۱06۵ 60۵۵ ۱۱۱۵8۵۵ 60۱۱۱۵2 ۱۵/۵۵ ۵۳۵۱۵1۵۲ 561511116 
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۲۵۱۱۵۱۵ ۵ ۱۹۵0۲ 2۲۵-۵65۱ 811 ۰ 


کون ملمنامط نو رهم1۲ ۱۷۵0۵1 15 ,عصمتماموظ ممتامز۱60ظ رومامصهنوظ ممتما۱۷ فصامن عج0تا ۲۷۷۲۵۲۵8 رمک 
و( 


* تاریخ دریافت مقاله ۹۵/۱۱/۲۳ و تاریخ پذیرش آن ۹۷/۷/۱ می‌باشد. 
(۱) نویسنده مسئول: استادیار گروه مهندسی عمران, دانشکده مهندسی. دانشگاه شاهد. تهران. 

90۵860:6۲ ۵ 2072 :اق۲ 
(۲) دانشجوی دکتری. دانشگاه آزاد اسلامی» تهران مرکزی. 


۹۸ استتتاج معادلات جدید برای تخحمین بيشینه شتاب و سرعت زمین ناشی از زلزله 


مقدمه 
زمین‌لرزه به‌عنوان یک رویداد غیر قابل پیش‌بینی» همواره 
بشر را در معرض خطرات و صدمات جانی و مالی قرار 
داده‌است. ازاین‌رو خطرات و عواقب ناشی از زمین‌لرزه 
ناگزیر اما کاهش اثرات منفی آن امکان‌پذیر است. که این 
مهم نیازمند برآورد خطر زمین‌لرزه به‌ویژه در نواحی 
لرزه‌خیز است. ازجمله مهم‌ترین بخش‌های مطالعات 
ارزیابی خطر زمین‌لرزه برای یک ساخت‌گاه. پیش‌بینی 
پارامترهای جنبش نیرومند زمین است که به‌وسیلة روابط 
موسوم به روابط کاهندگی يا روابط پیش‌بینی حرکات 


ممتامن۳۵0 ممتاما .. 0صتا06۲۵)) 


نیرومند. زمین 
(معمتاتو۳ صورت می‌گیرد که گزینش مناسب‌ترین و 
کاراترین روابط کاهندگی» یکی از مهم‌ترین اهداف 
محققان در تحلیل خطر لرزه‌ای یک منطقه به‌عنوان یکی 
از اصلی‌ترین اجزای تحلیل خطر لرزه‌ای و یکی از مراحل 
بنیادی در طراحی سازه‌های مقاوم دربرابر زمین‌لرزه است. 
پارامترهای مهندسی یک زمین‌لرزه را می‌توان به دو رده 
اصلی پارامترهای مربوط به حوزه زمان مانند بیشینه شتاب 
زمین (۸66۵۱62000 ۲0000 ۳۵۵) و بیشینه سرعت 
زمین (۷۵۱۵00 0۲0000 [۳۵۵) و نیز پارامترهای مربوط 
به حوزه پاسخ مانند شتاب شبه طیفی 5066721 9600ظ۳) 
(جمتاحه [ععع خر تقسیم نمود که می‌توان برای ارزیابی 
ریسک سازه‌ها از آنها بهره برد. با وجود اين که پارامترهای 
طیفی کارامدتر از پارامترهای حوزه زمان هستند [1] 
به‌دلیل مستقل بودن پارامترهای حوزه زمان به‌کارگیری 
آنها آسان‌تر است؛ بنابراین در پژوهش‌های خطر لرزه‌ای 
0۵ و ۲0۷ بیشتر مورداستفاده قرارگرفته‌اند. به‌منظور 
برآورد این عناصر می‌توان از روش‌های مختلفی مانند 
بررسی درجا و يا مدل‌سازی فیزیکی استفاده نمود که 
فخهو بستبار کشتوده تران فزسا ق پرش یه ات21 زا 
معادلات پیش‌بینی حرکات زمین (ع01۳۳) نقش 


برجسته‌ای در آنالیز خطر لرزه‌ای بازی می‌کنند. ۳0۸ 


نشریه مهندسی عمران فردوسی 


و۳0۷ پارامترهای حوزه زمان حرکات زمین هستند که 
به‌ طور متداول مورداستفاده قرارگرفته‌اند. شاخحص‌ترین 
جنبه‌های لرزه‌شناسی که پارامترهای حرکات زمین را 
تحت تأثیر قرار می‌دهند شامل اثر منبع» اثر مسیر و اثر 
سایت است. اثر منبع به پارامترهای گوناگونی هم‌چون 
سطح افت تنش در رویداد زمین‌لرزه. مکانیسم گسلش و 
جهت گسلش مربوط است که اثر اين پارامترها یکسان 
یت ار ال بات موروطر از کی یس زاره 
و برای آن در ادبیات فنی موضوع تعاریف مختلفی (مانند 
یک وی ال یل گیگ امه تور 
و 
به‌صورت کلی درنظر گرفته می‌شود (به‌عنوان مثالء خاک 
یا سنگ) و در برخحی دیگر از سرعت موج پرشی در خاک 
په‌عنوان شاحص آثر سایت استفاده می‌کنند.به‌طور کلی» 
5 نوین عموماً ازنظر بزرگی زمین‌لرزه» فاصلة منبع 
تا سایت. شرایط ژئوتکنیکی سایت و مکانیسم گسلش 
بررسی می‌شوند [13-5 اما پارامترهای فیزیکی موثر دیگر 
مانتل اف تنم اتشار کسککی: راسالیری: و رشاو 
غیرخطی خاک به‌دلیل پیچیدگی بسیار زیاد در توسعة 
معادلات پیش‌بیتی مورد استفاده قرار نمی گیرند [67]. 

فا لین این بر وهسن وشن اتی بی صبته تایب 
زمین (۳0۸) و بيشینه سرعت زمین (۳0۷) با بهره‌گیری 
از الگوریتم درختی مبتنی بر رگرسیون موسوم به 15 
است که پس از تهية مدل‌های نهایی برای هر یک از 
انواع سل به‌منظور اعتبار سنجی, مدل‌های ارائه شده با 
۳ مدل بیشتر شناخته‌شده در این زمینه ( 200 ۳۴00۲6 
ممعصتلاه [8],منصع/:020ظ همه لاهطامرصصعت [9] و 
21 )6 ,تصطمجع0. [10]) که دارای اشتراک در پارامترهای 
مورد استفاده در مدل‌سازی بودند مقایسه گردید و برای 
تحلیل عملکرد و اعتبان آنالیز پارامتریک و حساسیت 
برروی مدل انجام شد. 


سال سی و یکم» شمارُ چهار, ۱۳۹۷ 


علی درحشانی- علی صابری 


الگوریتم 15 
با پیشرفت در روش‌های داده‌کاوی به‌کمک رایانه. 
معادلات جدید برای پیش‌بینی حرکت زمین برپاية 
الگوریتم‌های ژنتیکی و ترکیبی ژنتیکی [10] شبکه‌های 
عصبی مصنوعی (۸(۷) [10,11] و دیگر الگوریتم‌های 
هزات فص یی [ و 12| اراقه تنل انست و اه ییا 
شبکه‌های عصبی مصنوعی ( لفتسعلا لمامتگنا۸ 
به‌طور گسترده‌ای برای پیش‌بینی معادلات 
بعرکته خی نان انیت ابا ای رون اه 
بسیاری از دیگر انواع روش‌های الگوریتم‌های هوش 
مصنوعی (معممونااعاص1 امژهتلتا۸) شفاف و قابل‌فهم 
نیست. ازاین‌رو به‌عنوان یکی از اهداف اصلی این 
پژوهش. به‌منظور یافتن روابط ساده‌تر از روش رگرسیون 
وه وا یاه لآشتی انم نی نکن اه 
موفق‌ترین و رایج‌ترین روش‌های مدل‌سازی‌های 
داده‌محور است که برای نخستین بار توسط حداجند0 
[14] به‌صورت مدل درخت تصمیم موسوم به ۷15 برای 
پیش‌بینی داده‌های پیوسته ارائه شد و ۷۷۵0۵ 200 ۱۷۷1۱660 
[15] این مدل را توسعه دادند. برخلاف مدل‌های درحت 
تصمیم معمول که کلاس يا رده‌های گسسته را به‌عنوان 
مخروجی ارائه می‌کنندء مدل 1۷15 مدل خطی چندمتغیره را 
راع د اما وراه کرو تال شرشی مس ساره بای 
ساختار مدل‌های درعت تصمیم‌گیری دارای ۳ مرحلة 
صلی است؛ ایجاد درخت. هرس کردن آن و ساده‌سازی 
درخت. در گام ساختن درخت. از یک الگوریتم 
ستنتاجی يا معیار تقسیم (انشعاب) برای ساخعت یک 
درخت تصمیم استفاده می‌شود. معیار بخش کردن برای 
لگوریتم مدل 145 ارزیابی انحراف معیار مقادیر کلاسی 
ست که به‌عنوان کمیتی از خطا به یک گره می‌رسد و 


کاهش مورد انتظار در این خطا را به‌عنوان نتیجهٌ آزمون 


هر پارامتر یا صفت در آن گره محاسبه می‌کند [16]. در 
شکل (۱) نمونة جداسازی درخت نشان داده شده‌است. 
کاهش انحراف معیار (]8) از رابطةّ (۱) به‌د.ست 


سال سی و یکم» شمارة چهار, ۱۳۹۷ 


۹۹ 


59۲ - 2009 - 400 ۷۳ 


ت ۷۵5 
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شکل ۱ جداسازی فضای ورودی و پیش‌بینی توسط مدل درختی 


برای داده‌های جدید. 2) جداسازی فضای ورودی ۶2 1 به‌وسیلةً 
الگوریتم درختی 5 و «) پیش‌بینی مقدار جدید توسط مدل 


درختی [17] 


که در آن ۲ بیانگر یک سری نمونه‌هایی است که به 
ی وا ای ۱ تیان کر مموتههاتی ات کر لسن 
فروجین سری. پتانسیلی. .وا دارند ق 50بیانگر انضرات 
معیار است. به‌دلیل فرایند انشعاب. داده‌های قرارگرفته در 
گره‌های فرزند. انحراف معیار کمتری نسبت به گره مادر 
داو ند و ینار آنر تالف نز هد مش اتسار همه 
انشعاب‌های ممکن» 2۷5 صفتی را گزینش می‌کند که 
کاهش مورد انتظار را بیشینه کند. این تقسیم بیشتر ساختار 
شبه‌درختی بزرگی را تشکیل می‌دهد که سبب بیش‌برازش 
(0۷6۲-۳10108) می‌شود. برای غلبه بر مسئلة بیش‌برازش: 
درخت تشکیل‌شده باید هرس شود. این کار با جایگزینی 
یک درخت فرعی با یک برگ انجام می‌شود. بنابراین؛ 


گام دوم در طراحی مدل درختی شامل هرس کردن 
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۱۰۰ استتتاج معادلات جدید برای ت- نحمین بیشینه شتاب و سرعت زمین ناشی از زلزله 


درخحت رشدیافته و جایگزینی درختان فرعی با توابع 
رگرسیونی حطی است. این روش ساخت مدل درختی. 
فضای پارامترهای ورودی را به نواحی يا زیر فضاهای 
کوچک‌تر تقسیم می‌کند و در هرکدام از آنهاء یک مدل 
رگرسیونی خطی برازش می‌دهد. بعداز اين‌که مدل خطی 
به‌دست آمد, برای کمینه کردن خحطای تخمین با حذف 
کردن پارامترها؛ ساده‌سازی مدل انجام می‌شود. در مدل 
5 از یک جستجوی حریصانه برای حذف متغیرهایی 
که مشارکت کمی در مدل دارند» استفاده می‌شود [14]. 


پارامترهای ورودی 
در این پژورهش بر پایة مرور ادببات فنی» 3 پارامتر در 
مدل‌سازی معادلات پش‌بینی ۳0۸ و ۳0۷ به کار 


رفته‌است. فرمول‌بندی درنظر گرفته‌شده بدین‌صورت 


۳0۸۵ 


)۲( ( ۲ مووا زا ٩)‏ رومم 


که در آن ۷1۷ بزرگی زمین‌لرزه» «یع1 نزد یک‌ترین 
فا صله افقی تا تصویر عمودی صفحد گسیختگی گسل 
(فا صلة جوینر- بور)؛ ۷:0 متو سط سرعت موجبر شی 
در ۳۰ متر بالایی سایت و ۳ مکانیسم گسلش شامل 
گسل‌های نرمال, امتدادلغزه معکوس و معکوس مایل 
است که باتوجه به زاوية ریک و مطابق جدول (۱) و 


کل (۲) تفکیک می‌شود: 


پایگاه داده‌ها 
منبع داده های حر کات نیروم‌ند زمین به کاررفته برای 
تو سعه 011۳125 پایگاه داده‌هایی بوده‌ا ست که در پروژه 
عصنمهمصتعص ملقتمطه ۴‏ مرمع ۳۳۴۲۱۲۲۵ 
مامت ممتهمج0۵) هلا عاصمی ‏ دوعوم 
(متطاعمم5612)1 به‌و سبلهٌ 2 اه ,۳0۷۵ [18] جمع اووم 


شده بود. 
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جدول ۱ تقسیم‌بندی مکانیسم‌های کشاسن برپایةٌ زاوية ریک 
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شکل ۲ تقسیم‌بندی مکانیسم‌های گسلش برپایه زاویةٌ ریک 


پایگاه داده (7.3 ۷ ۳12816 0۸ دربرگيرندة 
داده‌های ثبت‌شده زمین‌لرزه‌های پوسته‌ای کم‌عمق در 
مناطق تکتونیکی فعال جهان است. پایگاه داده‌ها طیف 
کسترهه‌ای. از نززکی. و فاضله را پوشسشی می‌دهد, 
پژوهشگران بسته به نوع مطالعه ممکن است که از همه 
پایگاه داده‌ها استفاده کنند و پا آنالیزهایشان را به 
زیرمجموعه‌های گزینشی محدود کنند. در اين پژوهش. 
بخشی از پایگاه داده از آنالیز حذف شده‌است. مجموعةٌ 
داده‌های فاقد اطلاعات پارامترهای موردنیاز و هم‌چنین 
گاشت‌های تکراری موارد حذف‌شده از آنالیز هستند. 
درنهایت. از مجموع ۱ نگاشت. ۲۷۷۷ نگاشت برای 
انواع گسل‌های ذکرشده شامل گسل امتدادلغزه نرمال 
معکوس و معکوس مایل برای ایجاد مدل استفاده شد. 
متغیرهای پیش‌بینی کنندهٌ موجود در آنالیز »21 هنک م۷ 


و ۴ بودند و ۳0۸ و ۳0۷ پارامترهای حرکات زمین 


سال سی و یکم» شمارُ چهار ۱۳۹۷ 


علی درخشانی- علی صابری 


بودند که فرموله شدند. شکل (۳ توزیم داده‌های 
مور خانتفاهه رای با شک ال ماد لک پیت بش را تشاق 
می‌دهد. 

محدوده‌ها و آمارهای پارامترهای درگیر در فرایند 
مدل‌سازی در جدول (۲) نشان داده‌شده‌است. همان‌طور 


۳ 


هستنك. کتت ف بودن محدوده‌ها این امکان را می‌دهد تا 
بتوان در مسائل مختلف و با اطمینان بیشتری از معادلات 
پیش‌بینی نهدستت أ سم استفاده نمود؛ لذا این کسفره کی 


نقش مهمی در آنالیز خطر لرزه‌ای دارد. 


برای تجسم بهتر, داده‌ها به وسیلة نمودارهای 


هیستوگرام نشان داده‌شده‌اند. 


که در این جدول دیده می‌شود دامنه‌های متغیرهای 


پیش‌بینی کننده» به‌ویژه برای بزرگی و فاصله نسبتاً گسترده 


۰۹2 ۳ 3 
۵ دی ۵ 2 ۵ ۵ ی نت یج ت 
وه ۱ ۳ * ویبوااا0 - ا,ملا . * 
وق و4 هن ۵ 5 ۲۳ 
عباواا۵0 - ۲٩6۷6۲۵‏ 
ما٩‏ - ٩۱۲۱6‏ 
۳ ۲ ۲ 1 5.0 
1000 100 10 1 01 001 


دز 
شکل ۳ توزیع داده‌های مورداستفاده برای توسعه معادلات پیش‌بینی تفکیک شده برپاية نوع گسل 


جدول ۲ آمار توصیفی از متغیرهای مورداستفاده برای ساخت مدل 


(628/560) ۳۵۷ | (ع) ۳0۸۵ | (5ل) وجی ۷ | (ص) وک ]۷ ۳2۵۵۵۵۲ 
938 0068 2900 334 6355 صهع]۷( 
024 0.00 3228 096 0.01 ۶۲ مها 
19 004 25.42 6349 6.30 صهنع]۷ 
٩1200270 06۷121008 9 305 1129 011 1265‏ 

۱/۹ 0.35 | 2555.97 | 3 001 1600214 

۱۹۹ -068 219 1022 2392 1149 

328 143 212 138 03 ۱۱۹۵ 
11104 16 19978 1 2.70 م۴ 
010 0.00 11635 0053 5.20 1( 
111314 16 2016-13 4 ۲ 7.90 ۱/۹ 


سال سی و یکم شمارة چهار» ۱۳۹۷ نشريه مهندسی عمران فردوسی 


۱۳۲ استنتاج معادلات جدید برای تخمین بیشینه شتاب و سرعت زمین ناشی از زلزله 


10096 ۴ 4 ۲ تِ 18 1032 1200 10096 1400 
9090۰ 9034 
20 1000 20 1200 
7096 200 7096 100 
ترات ات 800 
5096 500 5094 
4096 4096 600 
2096 ۵« موم 400 
2096 2090 
10896 1 1 بت كِ 109 200 
09۶۰ «___ ستت 0 09 0 
۶ 350 250 ۰ 300 ۰ 200 0 50 ۰۶ 750 700 650 600 . 550 عععا 
() دز ۱۷۳ 
1423 0 ۰ 10096 1600 
100 _ 100 
0 پوو2 1200 
1000 5096 1000 
۳9۹ مومع . 5096 800 
4096 
500 وت 600 
400 2096 190 ٍ 
6 200 109۰ 35 200 
د ۳ - 0 8 8 # . 
مااک5:۳۷ 1850 1600 1350 1100 850 600 ۰ 350 عععا 
ععها6 ۱۷۶۵۵06 (۱۱/5) ۷530 
600 
500 
400 
300 
200 
100 
0 
(0۳/5) 26۷ 
46 یتست ۲۲0 ۳۳ 
شکل ۶ هیستوگرام پارامترهای مورداستفاده در مدل‌سازی 
پیش پردازش داده‌ها می‌گیرد؛ درذها یت ع1 از هر دو سوی رابطه حذف 
خروو جی الگوریتم 5 به صورت جمع و تفریق می‌گردد. مراحل مذ کور به ترتیب در روابط (۶ (۵) و 
پارامترهای ضریب‌دار و به‌شکل زیر انتت: )1( نمایش داده شده‌است. 
20۸۵ 
)۳( و۲۷ مر ۵۶ ۳۳ < برجم ط بط +6 + وب - 1۱۳۵۸ 
۱ ۴ 5 ۱۳۳0۷ 
که در آن ۰ 8 و ۷ ضرایب پارامترها و ) مقدار 6 + وچی۷ 10 * ۷ 
در این پژوهش برای آسان‌تر شدن فهم معادله و 
اه که پارامتر چه تاه اه ها ۳۸ ۲6۸۱ 18 
و ۳0۷ دارد. نخسست از همه پارامتر های ورودی و ) 6 + و۷ 10 * ۲ + دز 13 + ۵ 


نشریة مهندسی غمراآن فردوسی سال سی و یکم شمارة چهار؛ ۱۳۹۷ 


علی درحشانی- علی صابری 


20۸۵ 
20۷ 


0( 6 2 + ویو + تورنط * ترا 

پس از محاسبة 10 همه پارامترهای موردنظر ورودی 
و خحروجی ۲۷۷۷ داد گزینش‌شده در نخستین گام از 
پیش‌پردازش. داده‌ها برپایة نوع گسل از یکدیگر تفکیک 
کردیه یرای یقن هوریگه ان فضیوت داهان 
تفکیک‌شده پرحسب نوع گسل, ۸۰ درصد داده‌ها به‌طور 
تصادفی به زیرمجموع آموزشی و ۲۰ درصد داده‌ها به 
زیرمجموعة آزمایشی احتصاص یافت که در جدول (۳) 
شمار داده‌های اختصاص‌یافته برای هر یک از این 
زیرمجموعه‌ها نشان داده شده‌است. داده‌های یادگیری 
برای آموزش الگوریتم آماده شد و به‌وسيلة آنها مدل 
ساخته شد. داده‌های اعتبارسنجی وی مد شودست امه 
بااستفاده از داده‌های آموزشی قرار گرفت و از آن برای 
مشخص کردن توانایی تعمیم مدل‌ها استفاده شد؛ 
به‌اين ترتیب. هم داده‌های یادگیری و هم داده‌های 
اعتبارسنجی در فرایند مدل‌سازی. لحاظ شدند. برای 
ارزیابی کارکرد مدل‌های به‌دست‌آمده به‌وسیل 0۷15 
داده‌های آزمون که هیچ نقشی در ساخت مدل‌ها نداشتند. 
در مدل نهایی قرار گرفتند و مقادیر خطا و ضریب 
همبستگی که درادامه توضیح داد شله‌است:به‌ذست امد: 
به‌منظور به‌دست آوردن یک دسته‌بندی مناسب از داده‌هاء 
چندین ترکیب از مجموعه‌های آموزشی و آزمون مدنظر 
قرار گرفت. گزینش به‌گونه‌ای بود که بیشینه. کمینه. 
میانگین و انحراف استاندارد پارامترها. در مجموعه 


داده‌های آموزشی و ازمون سازگار باشد. 


جدول ۳ تعداد داده‌های آموزشی و آزمایشی در هر یک 


از انوا اع کسل 
1۳۵ | ل101۸ دعصشت ۷۲۳۵۲۱۵۸۵۲۲۹۲[ 
130 519 649 طلاو - ]۴ 5 
13 53 6 ۱9۴ 
284 1139 1123 ۷۳۵ ۴۳۴ 
125 511 639 0۳۴1/10/۴ - 1۱۳۴۱۷۴۳۴۵۴ 


سال سی و یکم» شمارة چهار, ۱۳۹۷ 


اعتبار سنجی 
برپایة مرور ادبیات فنی موضوع سه پارامتر عمده که در 
اکثر مطالعات به‌عنوان پارامترهای اعتبار سنجی مدل‌ها از 
آنها | ستفاده شدها ست. عبارتند از: ضریب همبستگی 
(اصمنء‌تلگهه) «منلقآم0۳)): میانگین ريشتء مربعات 
خطا (02ظ متهیاو5 مهع/1۷ ا0) و میانگین عطای 
مطلق (0۲ عادام‌دطاظ «2ع۷) که به صورت زیر 


تایه می تور ند: 


تسنیا فرص یط 


1 2 
86۳ < 2 ۳ 


1 
۸( ان - رط ,<< ۱۸۸۲ 


0 (0 


که در آنها زنط و زا به‌ترتیب. مقادیر خروجی واقعی 
و پیش‌بینی‌شده برای خروجی ام و زو :آ به‌ترتیب» 
میانگین مقادیر خروجی واقعی و میانگین مقادیر 
ی تن که دز تشوتدها و کیزاه تیرتها اس 

دقادیر, بالاتر 06و سقادیر تایی تر 1452و 
۲ تشان‌دهندة یک مدل دقیق‌تر است. 


ارائه و تحلیل مدل 

در این بخش مدل‌های به‌دسست‌آمده ب‌روش الگوریتم 
درختی 1۷15 برای بیشینه شتاب زمین و بیشینه سرعت 
زمین که برحسب نوع گسل تفکیک‌شده به‌همراه نمودار 
دوهی پیتو‌بیتی ] مدهآشسنسنت.. میسین نموداز, مقاذیر 
اندازه گیری‌شده دربرابر مقادیر به‌دست‌آمده به‌وسیلهٌ مدل 
برای مجموعه داده‌های آزمایشی و مجموعه داده‌های 
اعتبارسنجی رسم شده و با چند مدل دیگر که دارای 
اشستراک در پارامترهای ورودی مدل هستند مقایسه 
گردیده‌ا ست. سپس برای هر بخش تحلیل حساسیت و 
مطالعة پارامتری انجام شده‌است. 


نشريه مهندسی عمران فردوسی 


۳2 استتتاج معادلات جدید برای تخمین بیشینه شتاب و سرعت زمین ناشی از زلزله 


برای یا به کار گیری» روابط پیش‌بینی بیشینه نشان داده شده‌است و درادامه رابطه‌های ارائه‌شده 
شتاب زمین (۳0۸۵) و بيشینه سرعت زمین (۳0۷) در به‌صورت جدول (4) شامل توان‌های پارامترهای ورودی 


قالب نمودار درختی به‌ترتیب در شکل‌های (۵) و (") و ضرایب ثابت هر معادله آمده‌است. 


1 ۲۳/3100 6 > راز؟] 
0۵002 6 < از 


8 > راز 


3 > وز] 
8 << از 


۷ > 7 


3 < راز 
7 < ۸/۷ 


5 > از 


5 > «ز] 
5 < ور 


9 «متادنو] 5 < راز 
0 0 و22 > وز 
41 0۱31100 ] 9 < راز 


شکل ۵ چارت روابط پیش‌بینی بیشینه شتاب زمین ۳0۸ 


۴0۳/31003 


4 ۲0۲31100 
منک 5۱۲۱۷ 


۲90311005 


۴0۵31100 6 


۱ا۴ 


20۱۳03107 


۲0۵31100 8 


(۱-۵ ۹ 


عون6-00۱ 8۵۷۵۲6 


5۲۱۷۵ 0 


دتوترت را 


۴۶0۱۷31100 4 


۴0۱۷۵1100 15 


۴3۱ ۷ 


5 > ز9 


۲0۱۵110۳7 6 
[| < 5 


۴۲0۱۷21100 17 


56۷6۲56- 
(3 


[| > 2 ۲0۱۷21107 18 


2 < و9 


۴0۱۷3۲100 9 


شکل 1 چارت روابط پیش‌بینی بيشینه سرعت زمین ۳0۷ 


نشریه مهندسی عمران فردوسی سال سی و یکم شمارة چهار» ۱۳۹۷ 


علی درحشانی- علی صابری 


مقایسه با مدل‌های دیگر 


ی ۵ ۷۵0 | ۳0۷۷۵۲ وک[ 
1.757۳-1 4 6- 
3.272۳2 91 5- 
1.653۳1 1 6- 
7.608۳5-2 1 027 
00-+]2.4441 7- 6 
1.666۳0 52.- 2.-- 
691۳-04 5 3 
2.1۳-5 1 4- 
1.022۳-1 7 0.- 
1.461۳1 3.- 2- 
1194۳700 4 5- 
441۳-1. [ 1- 6 
1.191۳2 4 69- 
143۳-04 0.042- 5 
2.475۳3 92- 0- 
02-+7.161۳1 0.80- 05- 
5145۳۷2 28- 0- 
2.02۳02 4 1- 
2.045۳3 017- 472.- 


برای اعتبارسنجی. مدل‌های ارائه شده برای بيشینه شتاب 
زمین (۳0۸) و بیشینه سرعت زمین (۳۷) برپایة 
ضریب همبستگی (6» میانگین ريش مربعات خطا 
(9۳() و میانگین خطای مطلق (0۷۵۳ با ۳ مدل 
ارائه شده به‌ وسیله ۸۶09 200 00۲6ظ 
٩070:8212 ]8[‏ صرح ا[مطامرصفت [۲9] و ل2 اه رتصطملصون. 
[10] که دارای اشتراک در پارامترهای استفاده‌شده در 
مدل‌سازی بودند. مقایسه شد. از برتری‌های مدل 
ارائه‌شده در این پژوهش می‌توان به محاسبه مستقیم 
۸ و ۳0۷ به جای 10 اين مقادین سادگی مدل و 
برتری از نظر ضریب همبستگی (06) در بیشتر موارد. 
میانگین ری شذ مربعات خطا (1808515) و ميانگین حطای 


سال سی و یکم» شمارة چهار, ۱۳۹۷ 


۵ ۷۱ | صمتافتاوظ | 00۲۳۲ 

1 2352 

2 32-138 

3 0444 

4 0444 

5 00 

۳0۸۵ 6 090 
7 2-2-93 

8 602 

9 2-34 

10 0529 

11 109 

12 6 

13 1-114 

14 692 

15 358 

۳0۷ 

16 0512 

117 0512 

18 4210 

19 0198 


مطلق (۷۵۳) در انواع کی« هم در مجموعه داده‌های 
آموزشی و هم در مجموعه داده‌های آزمایتتی اشاره 
نمود. مقایسه برروی داده‌های آموزشی و آزمایشی 
استتفاده شده در فرادند ابداع مدل های حاضر 
انجام‌شده‌اسست. همان‌طور که درادامه خواهد آمد 
ملاحظه می‌شسود که در قریب به‌اتفاق موارد. برتری با 
مدل به‌د ست‌آمده با روش 1۷15 است. هم‌چنین, اگر یک 
مدل 00<0.8 ارائه دهد و مقادیر خحطا ( به‌عنوان‌مثال» 
۲ و ۷۸۲) حداقل پاشتله یک همیست کی فقوت 
میان مقادیر پیش‌بینی شده و اندازه‌گیری شده وجود دارد 
زا من مان فان بسن شرف ور 
گرفته شود. مدل های ۷5 به‌دست‌آمده با 00 بالا و 
مقادیر 1۷/5 و ۷1۸ کم به‌میزان قابل‌ملاحظه‌ای قادر 


نشريه مهندسی عمران فردوسی 


۱۰۹ استتتاج معادلات جدید برای تحمین بيشینه شتاب و سرعت زمین ناشی از زلزله 


به پیش‌بینی مقادیر هدف هستند [20]. () مقایسه شد. مشاهده می‌شود که در اکثر موارد 
مدل به‌دست‌آمده برای بیشینه شتاب زمین و بیشینه برتری با مدل 1۷۲5 است. 


سرعت زمین با ۳ مدل دیگر مطابق جدول (۵) و جدول 


جدول ۵ مقایسةٌ مدل ۷5 با ۳ مدل دیگر ۳0۸ برای داده‌های آزمایشی در هر یک از انواع گسل 


۳ 9 ٍ ۹1 ب ك_ ین 
2008 20017 اسب نذا 

0933 99617 0 39904 0 0(289988) ن 
00340 994 ۰ ۰320 0 00۳23332 ۱/۳2۹9 آ 
00567 (وز: ۱۳۱ 89 9 003150 ۴(۳۷۹۳ 
9106 00 ۰7 0 09505 0۳39905 00 
0364 00 09۰3130 0017 032163 1۵۳ 5 
0020956 ۰.9353« 0113 0038 ۱/۶ 
31 (۱0۳۳/۱ 8 0۳3۵317 09 
03913« 161 0 985« 0۳29121 ۱/۳2۹9 4 
00217 ۰055« 1-1346 9313 00 ۱/۶ 
239041 0 015 00 05142 0026910 9 
03927« 026532 165 00 23997( ۱/۳۹۹ ۹۵ 
00662 ۰99017« 00289016 003662 ۴(۹۳ 


جدول 7 مقایسةٌ مدل ۷5 با ۳ مدل دیگر 20۷ برای داده‌های آزمایشی در هر یک از انواع گسل 


0 
20 20017 200 

9 9 9015 ۱ 9914 0991 002991 

۹-1075 023136 0023 ا0۳۳ ۱/۳2۹9 آ 
1-640 2224 0۳ 394 0 915« ۴۳۷/۹۳ 

۵ 0 22651 00 899244 0021922 «9939 

5 ۱/۳2۹9 00 3919 00 162 0 23346 2-924 
۴(۳/۹۳ 00 239906 0 2۰46 02500 69 

0 119 00 3915 019 0 

۴ ۱/۳2۹9 0 133 05059 023405 2-25 
۴۳۷/۹۳ 031448 0۳23690 0 26958 14 

۵ 0۳ 0۳2 003916 0289500 

۹۵( ۱/۳2۹9 026366 00 932 0۳ 96 61617 
۱/۵ 10323 02038 023990 99138 


| سال سی و یکم شمارة چهار, ۱۳۹۷ 


علی د رحشانی- علی صابری 


تحلیل حساسیت برای درک آن است که هر یک از 
پارامترهای ورودی ازجمله بزرگی زمین‌لرزه. فاصلهةٌ 
جوینر- بور و متوسط سرعت موج برشی در هر یک از 
انواع گسل, در ایجاد مدل چه اندازه تأثیرگذارند. به اين 
منظور به‌صورت جداگانه. برروی هر یک از مدل‌های 
پیش‌بینی بیشینه شتاب زمین ۲0۸ و بیشینه سرعت 
زمین ۳0۷ تحلیل حساسیت انجام شد. روش انجام 
تحلیل حساسیت برای هر پارامتر ورودی در یک گسل 
به این شکل است که از داده‌های آموز شی موردا ستفاده 
در ایجاد مدل. پارامتر ورودی موردنظر حذف شده و 
سپس با شرایطی دقیقاً همسان که برای ایجاد مدل اصلی 
ا ستفاده شده! ست. مدل سازی انجام می شود. پس‌ازآن 
داده‌های آزمایشی برروی مدل جدید پیاده می‌شسود و 
ضریب همبستگی (06) میانگین ريشة مربعات خطا 
(۷55) و میانگین خطای مطلق (۷۸۳۲) محاسبه 
می‌گردد. در جدول (۷) نتایج تحلیل سا سیت بي شینه 
شتاب زمین در هر یک از انواع گ سل برای پارامترهای 
۷۷ ۷۰:0 و «ن3 مطابق روشی که گفته شد ارائه 


شده‌است. مشاهده می‌شود که در همه انواع کته 
پارامتر فاص له «زک1 تأثیر زیادی برروی مدل دارد و در 
نبود اين پارامتر ورودی» ضریب همبستگی () 
میانگین ریِشْه مربعات خطا (101517) و میانگین خطای 
مطلق (1۸5) دسستخوش دگرگونی بس یار ز یادی 
شده‌است. در گسل‌های نرمال و امتدادلغزه پس از پارامتر 
فاصله و1 که بیشترین تأثیر را در ایجاد مذل دارد؛ 
پارامتر ۷1 به نسبت ۷,۵ تأثیر بیشتری در مدل دارد. در 
گسل‌های معکوس و معکوس مایل برخلاف گسل‌های 
نر مال و امتدادلغزه پارامتر ۷۰0 به نسسبت »1۷ تأثیر 
بیشتری در مدل پیش‌بینی بیشینه شتاب زمین دارد. 

تحلیل حساسیت بي شینه سرعت زمین نیز در هر 
یک از انواع کشت برای پارامتر های ب«]۷ ۷,۵0 و دنک 
محاسبه شده‌است. در جدول (۸) دیده می شود که مانند 
تحلیل حساسیت بیشینه شتاب زمین در هم انواع گسل» 
پارامتر فاصله «ز18 تأثیر بسیار زیادی برروی ایجاد مدل 
پیش‌بینی بیشینه سرعت زمین دارد. پس از پارامتر فاصلة 
«ی پارامترهای »6 و ۷,50 به‌ترتیب بیشترین تأثیر را در 
ایجاد مدل دارند. 


جدول ۷ تحلیل حساسیت بیشینه شتاب زمین در هر یک از انواع گسل 


۱/۵ ۱/۹۹ 09 
00567 00340 0933 
00۳06017 00312 08290 
0.0999 00371 00000 
0021 0۳0 009472 
956( 0 2036۹4«( 09106 
0.0909 0۳0 09903 
2159( 00۳102 01618 
00965 ۳۹06( 02 9059 0 
02617 00313 29531 
03962 0۳0۳ 0039545 
آزاز0۳۱ الز0۳0 39- 
0202 0.00320 029515 
0۰00662 0۳0 0029041 
00۳06817 0.0320 03899205 
0۳0 0040 03062 
020۳0663 2(314«() 009069 
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جدول ۸ تحلیل حساسیت بیشینه سرعت زمین در هر یک از انواع گسل 


۷۲( ۱/2۹۹ 0 
7.2040 65 09691 
93378 22204 08932 
1-2114" 66070 0.0000 
0-333 3.9900 904 0 
69 2-2-0924 00289939 
93373 2716 05975 
10236930 575 03993 
6324 2-547 929 0 
214 2-5 00289821 
3948 2-2111 2914 0۳ 
95527 2-2534 234(" 
۰-9157 2-26 002985 
9358 61617 08560 
از۱۳۳ 13047 0029585 
3614 1 11-3965 2337( 
9.9895 653 0289108 


بررسی پارامتریک 
در این پژوهش به‌منظور بررسی توانم‌ندی معادلات 
پیش‌بینی آنالیز پارامتریک برروی داده‌ها و مدل ارائه‌شده 
انجام شد. روش کار به این صورت بود که برای هر یک 
از مدل‌های پیش‌بینی بیشینه شتاب زمین ۳0۸و بیشینه 
سرعت زمین ۳0۷ تنها یک متغیر پیش‌بین ی کننده تغییر 
می‌کند در حالی که دیگر متغیر های ورودی در مقادیر 
میانگین مجمو عه داده های ارموق تاشت که" داشتقه 
می شوند. بدین صورت با نهادن این متغیرها در معادلات 
پیش‌بینی 4فسلعت | سل بهز یله قل: و۷1 و متسه 
مقادیر بیشینه شتاب زمین 0۸و بیشینه سرعت 
زمین 70۷ می توان نمودارهای برازش این مقادیر را 
دربرابر پارامتری که ثابت نبوده‌اسست. رسم نمود. 
تواماش بالات اراتزه امتح اد که شاوی 
پیش‌بینی‌شده چه اندازه با رفتار فیزیکی حاکم بر سیستم 
موودبروسی عم‌تغوانی ذازد تفیین شنه‌است 211 شکل 
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۶ 6966وهه طز 6ع0 ۷۲۵061 وعفلن صعنصهطهع]۷ 
1۷1۷ 
۱۹۱ 
۹ 
۷۹30 
1۷1۷ 
۱۱۹ 
۹ 
۷۹30 
1۷1۷ 
۱۵2۵ 
۹ 
۷۹30 
1 - 16۷76۲56 
« عنوناطاه 
۷۹30 


(۷) روند پیش بینی‌های مدل ۳0۸ را با تغییرات 
پارامترهای »2۷4 1 و ۷۰:0 نشان می‌دهند. نتایج آنالیز 
پارامتریک نشان می‌دهد که ۳0۸ همواره با افزایش ۷۲۷ 
افزایش می یابد و با افزایش «زک و ۷,50 کاهش می‌یابد 
که نان می د هد ذتایج آ نالیز پارامتر یک از د ید 
زمین‌لرزه‌شناسی همگی مورد انتظار بوده‌اند [22 ,3-5]. 

نتایج تأیید می‌کند که معادلات پیشنهادی برای 
بی شینه شتاب زمین ۳0۸ توانمند ه ستند و می‌توانند با 
اطمینان به‌منظور پیش گویی در مطااعات خطر لرزه‌ای 
مورداستفاده قرار گیرند. 

ش کل (۸) گرايش مدل پیش‌بینی ۳0۷ را به 
تغییرات پارامترهای لرزه‌ای» ۷ از و ۷۹30 نشان 
می‌دهند. نتایج آنالیز پارامتریک نشان می‌دهد که ۳0۷ 
همواره با افزایش ۷1۷ افزایش می‌پابد و با افزايش 18 
و ۷۰30 کاهش می‌یابد که نشان می‌دهد نتایج آنالیز 
پارامتریک با دیدگاه‌های زمین‌لرزه شناسی سازگار است 
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علی درخشانی- علی صابری ۱۰۹ 


و تائید می‌کند که معادلات پیشنهادی برای بیشسینه به منظور پیش گویی در مطالعات خطر لرزه‌ای 
سرعت زمین ۷ قوی هستند و می‌توانند با اطمینان مورداستفاده قرار گيرند. 
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شکل ۷ آنالیز پارامتریک برای بیشینه شتاب زمین برپایة مدل 5 
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شکل ۸ آنالیز پارامتریک برای بیشینه سرعت زمین برپای مدل 15 
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۱۱۰ استتتاج معادلات جدید برای ت- تتحمین بٌ بيشینه شتتاب و سرعت زمین ناشی از زلزله 


نتیجه گیری 

ها رن مها سس ها شتا 
زمین لرزه و به‌منظور پیش بینی پارامتر های حر کات 
نیرومند زمین ازجمله بیشینه شستاب زمین (۳0۸) و 
پیشینه سرعت زمین (۳0۷) به‌عنوان یکی از مهم‌ترین 
مراحل مطالعات ارزیابی خطر زمین لرزه. از یک روش 
داده‌کاوی مبتنی بر رگرسیون درختی موسوم به 115 بهره 
برده شد. 01۷1۳۲۳۵ پیشنهادی براساس و ۳/۳۷ 
نگاشت زمین‌لرزه از پایگاه داده‌های گسترده‌ای شامل 
۱ نگاشت زمین‌لرزه از سراسر جهان, جمع آوری‌شده 
در پروژ؛ ۳۷۳-0۸ ساخته شد. پارامترهای ورودی 
۷ بزرگی زمین‌لرزه «:*1 نزدیک‌ترین فاصلهٌ افقی تا 
تصویر عمودی صفحة گسیختگی گسل (فاصلةٌ جوینر- 
بور)؛ ۷:0 متو سط سرعت موج‌بر شی در ۳۰ متر بالایی 
سایت و 7 مکانیسم گساش به‌عنوان پارامترهایی کارامد 
برپایة ادبیات فنی موضوع در تهیة مدل‌های پیش‌بینی 
حر کات نیروه‌ند زمین به کار گرفته شد. برای تهیية 
مدل‌های 01۷۳105 داده‌های انتخابی براساس انواع ار 
تیا کنیا رها کصل دنه گس ود موی 
مایل که برپایهٌ زاویهٌ ریک تعیین می‌شود جدا گردید و 
داده‌های هر یک از انواع سل به زیرمجموعه‌های ۲۰ و 
۰ درصدی آموزشی و اعتبارسنجی تقسیم گردید و 
مدل با به‌کارگیری زیرمجموعه داده‌های آموزشی برای 
هر کش ایجاد شد. 01۷۳۳ پیشنهادی در این پژوهش 
برآوردهای قابل اعتمادی از مقادیر ۳0۸ و ۳0۷ ارائه 
می‌دهند و شرایط و معیارهای مختلف مد نظر برای 
اعتبار سنجی را ار ضا می‌کنند که از آن جمله می‌توان به 
موارد زیر اشاره نمود: 

معادلات پیشنهادی به‌طور چشم گیری ساده و 


فراجع 


جایگزین‌های کارامدی برای معادلات بسیار پیچیدة 


ارائه‌شده توسط دیگر محققان هستند. 


. نتایج آنالیز پارامتریک انجام‌شده برروی زیرمجموعه 


داده های آزمون با رف تار های موردانة ظار از د ید 
زمین لرزه‌شناسم هم‌خوانی دارند؛ به‌این‌ترتیب که 
۸ و۳0۷ همواره با افزایش ۷1۷ افزايش و با 


افزایش 89 و ۷530 کاهش می‌پابند. 


نشان می‌دهند که 0۷1۳129 استنتاج‌شله به‌نحو 
کارامدی رفتار فیزیکی | سا سی حاکم بر می ستم را 
درنظر می گیرند و تنها ترکیباتی از متغیرهای 
پیش‌بینی کننده ند سشند که به صورت تضادفی بهترین 
قانیت را با نتایج تجربی داشته باشند. 


. مدل های مربوط به ۳0۸ و ۳0۷ درمقایسه با چند 


مدل شناخته شده در ادیبات فنی موصوع که از لحاظ 
پارامترهای ورودی با مدل‌های پیشنهادی این پژوهش 


اششتواک-داشستت: ازنظر ضحر یب هفستکی (1 
میانگین ريشه مربعات حطا (۲(۷۲۲) و ميانگین 


خطای مطلق (۷۸۳) در بیشتر موارد برتر هستند. 


. نظر به این که 0۷1۳۳5 مبتنی بر 1۷5 برروی یک 


پایگاه داده‌های جامع با طیف گسترده‌ای از خواص 
توسعه یافته‌اند. می‌توانند با اطمینان برای اهداف 


طراحی به‌طور کاربردی مورد بهره‌برداری قرار گیرند. 


. افزون بر سسادگی روابط برخلاف اکثر مدل های 


پیش‌بینی حرکات نیرومند زمین محاسبةٌ بیشینه شتاب 
زمین (0۸) و بیشینه سرعت زمین (۳0۷) بدون 
استفاده از لگاریتم به‌عنوان یکی از برتری‌های مدل 


پیش‌بینی ارائه‌شده در این پژورهش است. 
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2011(۰) ,717-732 ۳۵۰ ,24 ۷۵۱۰ ,ع6ص۵عن00[1] 

عصتعده (۳۸) ممتاحمآهعع۸ 0۲۵۲0 »۲۵21 ۵ ممتامنلع۳۳. رتتمطتا2 رب مه بله۲ ویک رتقللزظ .6 رفمحصمط1 .و 
2013(۰) ,10 .۲ ر8 ,۷۵ رقلا0اا ۱۲60۳۷۵1 لقتناع۲ ۲))10121 خر 

م۱۷۱۵ مامت معط ما مهمتم۸ ۸02۷۵ مه رته‌طا۷۷۵ ,۷۷۰ 0۰ 200 رصه‌افظ بخ ات6216 ماه ۷ و.۳ 
,ما0۵ عمززممصزعص ر0مطام۷ ممتوععتعم عتامج0ه متتامجصم‌نهم- مهم 2 بط صمامامیظ ممتام۱۵01ظ 
2014(۰) ,1651-1668 .۳0 ,46 ,۲۷۵1 

0 هر عصلونا ومامصهنجظ ومتامص نام آهمتمصتوظ ۵۶ ممتام۱60۲ظ۰۳ رتصمل‌صعت .۲۲ مه 20 ۸27 .۲ خر 
۰ 51۳۱۵۱۱۲۵۹ > 0۱0۵۵۵۲ هه صصه 0ماهاناصل مه ۱۱۵۲۵0 تهتبهاظ لمتمتتانه عصتام‌دم ۱۷۲۵000 


89, ۴0۰ 2176-2194, )2011(۰ 


لمتمتکتانه وم ععصمعه‌کجمع کصئمز صفتلهتعنه اد صد ر ععععقهات فتامتهتاوم از مصتصعما رصقامتنای رب .ر . 


(1992) ,343-348 .30 رععصمع 10160111 


۲ ,۳۲0۵ طا ر 0128565 وتامتاصتاصمن) عصتام۲:601 ۲۵۲ قمع 1 ۷۲۵061 ۵۶ ممتامتا0ظ1. رصع ۷۷ ,۷۰ 220 ۷۷166۵ رب[ 
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سال سی و یکم» شمارة چهار, ۱۳۹۷ روتانس دوگ مرشر 


۱۱۲ استتتاج معادلات جدید برای ت- نحمین بیشینه شتتاب و سرعت زمین ناشی از زلزله 


1997(۰) ۲۵۵۲۱۱۱۸8۵ ۸/۵۵۱۵ 01 ۳۲۱۵۲۵۵۰ و۲۵۵ 

۶ ممام۲:۵0۲ ۶۵۲ مممتامتوظ و۱ اه ممتاه هن پتیج‌طام2 مها ریق مه ترمممصاعطاهظ رب بطه7 2 ریخ 
-1 .0( ر۱8 ۵۱ ۵۳۱۱۵۵۵ 0۴ 0۱۱۱۵ "رصطا01ع21 ۱۷5۲ مصاونا دمن ط10)مح- نام ۲۱00102 
2016(۰) ,21 

۵0۵1 0۶ ۰۸۵۵2010 رععزع۲ ریک یک صه رفصتامصمامی ر.ظ .1 رصقصصم‌امک ررظ رتاهع006ظ وی رال ین .۸ 
۵66۷01 هل موه ۷ رز متامت60:ظ 0۵ تحعاض ۲۵۲ عمامتصطمع1 عمتجم مصتطمه عطام 0صه وممع 1" 
2010(۰) ,262-274 0۰( و12 ۷۵۱۰ رکع0/۵۱۱۵ زن 0۱۱۵ ,1۵۲6۵۵ 0۶ 1۵000110 

ممتام]۱۷ مصنامتت) ۵۶ ممتاجهصون ۱۵ م۳۲ رععا9ام۲ بت 0ص2 مهدجه رل بتامنط ر.ظ ,۲۵0۷۷۵۲ .1۷۲ 
0۱۵۵۰۵۱۱۵۵ ۷۵۱۵۵ ۱و رقعصتمم۳۱۵ ظا ,0۷۵۲۹۵ ۸۵ ۳06۵ (هتاظ) و۱۷۵0 ومتاجیاهه‌ااض 
2۷6(۰) ۳و1 ۳۲۵1۵۵ 

6(۰) ,1۲08000 عصلامی :فمزهممصتعو 1 صا فمتافتاهای مه اتازطدطامظ رصان رل 0 

۲ ۶ ۳۲۵016108 ۲0۲ ۱۵061 ۵9۳۴ ۵۷۵1 ۵" بتتن ,.۲۷ 200 رم بیط بعظ۷۷2 ویک رعطقال بل بطق ۷۲ 
4 0۱۸۱۱۵ , "مصنطع۱۷]2 ما6۵ ۷ 9۱۵۵0۲۲ ون ۱860 ملصتامم‌صمن منصهع :0۵ 0۲ ماتصصنا زانازطاحجصها۲ 
2009(۰) ,962-969 ۳۴۵۰ ,168 ۷۵۱۰ ,عا111016۳16 

عصتنجءظ مطا عصتامنلع۳ 10۲ ۱۷۲۵061 0ععو-ظط۰۸ ره ععمک .۷۷ 20 رصتصصهرا بیق رهفهل ۷ بمتاک ۷ 
۰ ,36 ۷۵۱۰ ,و01661:1۱6ع) ۵۱۱ 0۳۱0۵/۵۲5 5011 1۷۵وعصم) ۱۲۵۱۵-12۵۵0 جم عصتام۲۵ متتگ ۵۲ مدومن 
(2009) ,503-516 

صد هامموک مقصممرعم؟ آماص۲۱۵۲۱20 ۵۲ ممتامنل۰۳۲۳۵ رتمصصحطمظ 1.۱ .ز فصع رصمومرحصو ررنا یک روعفهتصاصصخ ر لا ۷ 


1996(۰) ,371-400 .0 رد2 ۵1 ۷ ۱۳۱۵ 6۱۵۵۱۵ ک ع۱۵2۱۱8ع 1 ۵۳۱۱۵۱۵16 ۲۱۱۵۱۲0۵۵۵۲ 
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رنه هقی عی رال قرو رتور سال سی و یکم شمارة چهار, ۱۳۹۷ 


